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We theoretically investigate the equation of state and Tan’s contact of a non-degenerate three 
dimensional Bose gas near a broad Feshbach resonance, within the framework of large-expansion. 

Our results agree with the path-integral Monte Carlo simulations in the weak-coupling limit and 
recover the second-order virial expansion predictions at strong interactions and high temperatures. 

At resonance, we find that the chemical potential and energy are significantly enhanced by the 
strong repulsion, while the entropy does not change significantly. With increasing temperature, the 
two-body contact initially increases and then decreases like T~^ at large temperature, and therefore 
exhibits a peak structure at about ITco, where Tco is the Bose-Einstein condensation temperature 
of an ideal, non-interacting Bose gas. These results may be experimentally examined with a non¬ 
degenerate unitary Bose gas, where the three-body recombination rate is substantially reduced. In 
particular, the non-monotonic temperature dependence of the two-body contact could be inferred 
from the momentum distribution measurement. 

PACS numbers: 67.10.Ba, 67.85.-d 


I. INTRODUCTION 

Understanding strongly interacting Bos e gas es in three 
dimensions is a notoriously difficult quest [ 1 H 3 • Theoret¬ 
ical studies of these systems have been hindered by the 
absence of controllable theoretical approaches that can 
be used to describe their properties within certain errors. 
Although a formal field-theoretical description of weakly 
interacting Bose gases was developed more than half a 
century ag o by Lee, Huang and Yang and later 

by Beliaev based on the ground-breaking Bogoliubov 
theory [l^ . This theory is only applicable in the limit of a 
small interaction parameter, the so-called gas parameter 
nOg <C 1 - where n is the density and Ug > 0 is the s-wave 
scattering length - as a result of the perturbative expan¬ 
sion. When the gas parameter is extrapolated to infinity, 
each term appearing in the perturbative field-theoretical 
description diverges. To the best of our knowledge, a re¬ 
summation of these divergent terms remains unknown, 
even in an approximate manner. 

Experimental studies, on the other hand, have been 
hampered by atom losses from inelastic collisions. Un¬ 
like a strongly interacting Fermi gas, where the atom loss 
rate due to three-body recombinations into deeply bound 
diatomic molecules is greatly suppressed by the Pauli ex¬ 
clusion principle 0, at low temperatures an interacting 
Bose gas has a three-body loss rate proportional to of 
(i.e., the loss coefficient £3 ^ ha^jm, [isL 0), which 
grows dramatically when Cg is increased. Even in the ab¬ 
sence of inelastic collisions, for a strongly interacting Bose 
gas, the possibility of recombination into deeply bound 
Efimov trimers [l3| indicates that the system can be at 
best metastable. 

Due to these realistic problems, experimental studies of 
a strongly interacting atomic Bose gas near a broad Fes¬ 
hbach resonance have only been carried out very recently 
[T^-[0. The stability or lifetime of a unitary Bose gas 


with infinitely large scattering length was investigated 
with ^Li [0 and atoms [ 2 l| in the non-degenerate 
regime. It was found that there is a low-recombination 
regime at high temperatures and low densities, in which 
the loss coefficient saturates at L 3 ~ hX^^/m oc as 

predicted [0. Here, at high temperatures the thermal 
de Broglie wavelength Xds = [27:11^/(rnksT)]^/^ replaces 
the role of the s-wave scattering length Cg. The momen¬ 
tum distribution of a quantum-degenerate unitary Bose 
gas was also measured with ®®Rb atoms j^. These rapid 
experimental advances have trigged a number of inter- 
esti:M theoretical investigations on the unitary Bose gas 
[24 - K^ . focusing particularly on the universal Bertsch 
parameter the condensate fraction no at zero temper¬ 
ature and quenching dynamics. The predictions however 
are very different with each other, due to the absence of 
an efficient theoretical framework to handle the intrinsic 
strong correlations of a metastable unitary Bose gas. 

In this work, we aim to develop a non-perturbative, 
controllable theory of a strongly interacting Bose gas 
in its normal state, with an emphasis on the high- 
temperature low-recombination regime in which our the¬ 
oretical predictions might be efficiently tested in future 
experiments. Our description is built on an earlier inno¬ 
vative theoretical work by Li and Ho [0, who treated 
a repulsive Bose gas as a metastable upper branch (de¬ 
fined later) of an interacting Bose gas across a broad 
Feshbach resonance. By appropriately re-defining the up¬ 
per branch prescription |39| and using a non-perturbative 
large-7V expansion approach to remove the unphysical 
non-linear effect in pair fluctuations (dol - l^ . we overcome 
the large mechanically unstable area encountered earlier 
at low temperatures and therefore make Li and Ho’s 
idea practically useful at arbitrary temperatures in the 
normal state and arbitrary interaction strengths. Our 
improved theory is able to reproduce the path-integral 
Monte-Carlo results at weak couplin gs 14^ and the virial 
expansion at high temperatures |44l - l48| . In the strongly 
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interacting unitary limit, we calculate the equation of 
state and Tan’s two-body contact as a function of 
temperature. An interesting non-monotonic temperature 
dependence of the contact is predicted and is to be com¬ 
pared with future experimental measurements of the mo¬ 
mentum distribution. 

The rest of the paper is organized as follows. In 
the next section, we briefly introduce Li and Ho’s idea 
of the upper branch Bose gas and present the gener¬ 
alized Nozieres-Schmitt-Rink (NSR) method. The up¬ 
per branch is then appropriately defined through an in¬ 
medium phase shift. The large-expansion approach is 
adopted in order to overcome the unphysical strong pair 
fluctuations at large interaction strengths. In Sec. Ill, 
we first present the results for weakly interacting Bose 
gases and compare them with the available path-integral 
Monte Carlo simulations. We then discuss the equation 
of state and Tan’s two-body contact in the unitary limit. 
At sufficiently high temperatures, the results are com¬ 
pared with the virial expansion predictions. Finally, Sec. 
IV is devoted to the conclusions and outlooks. 


II. GENERALIZED NOZIERES-SCHMITT-RINK 
APPROACH 


A three-dimensional (3D) interacting Bose gas can be 
described by the imaginary-time action 


5 = 


/ 


drdx 





( 1 ) 


where tp^x), 'tp{x) are c-number fields representing the 
creation and annihilation operators of bosonic atoms of 
equal mass m at a space-time x = (x, t) . The imag¬ 
inary time T runs from 0 to the inverse temperature 
P = l/{kBT) and fi is the chemical potential. The inter¬ 
atomic contact interaction is parameterized by the bare 
strength Uq < 0, which has to be regularized by the two- 
particle s-wave scattering length Og via the relation, 

J_ = _- V — f2) 

Uo AirtPas V ^ 2ek’ 

k 


where V is the volume of the system and Ck = fi^k^/(2m) 
is the free-particle dispersion (i.e., kinetic energy). In 
experiment, the scattering length Og can be conveniently 
tuned by a magnetic field across a Feshbach resonance, 
to arbitrary values [li¬ 
lt should be noted that in our model action, Eq. m, 
the contact interaction is always attractive {Uq < 0), al¬ 
though the scattering length can change sign across the 
Feshbach resonance. This implies the pairing instabil¬ 
ity of two bosons and therefore the ground state of the 
system would be a mixture ofpairs and of the remain¬ 
ing unpaired bosonic atoms |51j . similar to what hap¬ 
pens for an interacting Fermi gas at the crossover from 
Bardeen-Cooper-Schrieffer (BCS) superfluids to Bose- 
Einstein condensates (BEC) |52 h^ . In the normal state. 


such a mixture can be described by using the seminal 
NSR approach 

Following the earlier work by Koetsier and co-workers 
1^ . we introduce a pairing field 

<(> (x, r) = DoV’(x, r) V'(x, r) (3) 

and decouple the interatomic interaction via the stan¬ 
dard Hubbard-Stratonovich transformation, with which 
the atomic fields appear quadratically and therefore can 
be formally integrated out. This leads to an effective ac¬ 
tion for the pairing field and, at the level of Gaussian 
pair fluctuations, results in the following grand thermo¬ 
dynamic potential: 

D = Do + (5D, 

Do = fcBT^ln(l-e-'5«‘‘), 

k 

sn = ksT ^ In [-F“^ (q,Wi)] , 

where ^k = Ck ~ The last equation is the contribution 
from pairs of bosons, which is characterized by the two- 
particle vertex function (or the effective Green function 
of pairs) with bosonic Matsubara frequencies 

vi = 2t:IT (Z = 0,±1,±2,---) [HI, 

IB {q, k) _^ 1 

~ 'fq/2+k ~ ^q/2-k 2ek_ 

( 7 ) 

Here nsix) = l/(e^“ — 1) is the Bose-Einstein distribu¬ 
tion function and the factor 7 b (q, k) = 1 + nB(Cq/ 2 -(.k) + 
riB('fq/ 2 -k) takes into account (in-medium) Bose en¬ 
hancement of pair fluctuations. By further converting 
the summation over Matsubara frequencies in Eq. ([6]) 
into an integral over real frequency and introducing an 
in-medium two-particle phase shift 

(5 (q,a;) =-Imln [-F“^ (q,a;-I-i0+)], (8) 

the contribution to thermodynamic potential from the 
bosonic pairs can be rewritten as 

1 r+oo 1 

= (9) 

To make the above integral meaningful, it is easy to see 
that, the phase shift at zero frequency w = 0 should 
vanish identically for any momentum q because of the 
Bose-Einstein distribution function. This is the so-called 
Thouless criterion, which is used to determine the onset 
of pairing superfluidity. 

We note that, within the NSR approach, the only pa¬ 
rameter in the imaginary-time action - the chemical po¬ 
tential /r - is to be determined by using the number equa¬ 
tion, 


F-i = 


m 


Airfflc 


-E 


( 4 ) 

( 5 ) 

( 6 ) 


1 d (Do -I- 5D) 
V ^ 


= no -I- Sn, 


( 10 ) 
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FIG. 1: (Color online), (a) The in-medium phase shifts for 
an attractive Bose gas S^ttiq = 0,^;) at the gas parameter 
nal = 1. The real and imaginary parts of the negative in¬ 
verse of the two-particle vertex function, —= 0, a;) are 
also shown, (b) The corresponding in-medium phase shifts 
5rep(5,(u) in the meta-stable upper branch at different mo¬ 
menta q — 0 (black solid line), kp (red dashed line) and 2kp 
(blue dot-dashed line). The phase shifts drep)? = 0, tu) at dif¬ 
ferent gas parameters na^ = 0.01 and na^ = — 1 are shown by 
gray cross and green stars, respectively. For all the plots, the 
temperature is fixed at T = 2Tco, where Tco — 0A36Tp is the 
condensation temperature of an ideal Bose gas. The chemical 
potential p is fixed to that of an ideal Bose gas at the same 
temperature, i.e., p = ii^^\2Tco) — —0.358£f. 

where n is the number density of the system, consisting 
of both the densities of atoms nq and of pairs 5n. 

For an attractive Bose gas near broad Feshbach res¬ 
onances, Eq. or Eq. oni) physically describes an 
ideal, non-interacting mixture of bosonic atoms with des¬ 
tiny no and pairs with density 5n > 0. With increasing 
strength of attractive interactions, the contribution from 
pairs, Eq. becomes more and more significant. As 
a result, the chemical potential decreases to the half of 
the binding energy, p —>■ —h'^/{2ma^)^ as required by the 
Thouless criterion 6 (q, w) = 0 (5ll |. 

A. In-medium phase shift for the ground state 

In Fig. UKa), we show the typical behavior of the in¬ 
verse vertex function F“^ and of the in-medium phase 


shift 5att for an attractive Bose gas with the gas parame¬ 
ter nal = 1 at T = 2Tco, where Tco ~ 0.436Tp is the con¬ 
densation temperature of an ideal Bose gas, measured in 
units of Fermi temperature Tp = ?i^(37r^n)^/^/ {2mkB) = 
Sp/ks- The phase shift jumps from zero to tt at the 
threshold frequency Wb(q), which signals the existence of 
bound states. Upon increasing the frequency beyond the 
scattering threshold, 

(q) = - 2m, (11) 

where the imaginary part of the vertex function be¬ 
comes nonzero, the phase shift decreases towards tt/2 
as w —> -|-c». Therefore, there are two contributions to 
the phase shift, originating from the bound states (at 
Wf,(q) < w < Ws(q)) and from the scattering states (i.e., 
uj > a;s(q)), respectively. It is clear from Fig. [T](a) that 
the phase shift, as an illustrated example, does not sat¬ 
isfy the constraint <5 (q, w = 0) = 0. This is because we 
have used an artificially large chemical potential, larger 
than the actual chemical potential, which has to be solved 
self-consistently by using the number equation (flUl) for an 
attractive Bose gas. 

B. In-medium phase shift for the upper branch 

It is interesting that although we are dealing with an 
attractive Bose gas, we may also obtain useful informa¬ 
tion about a repulsively interacting Bose gas, by treating 
it as a metastable upper branch of the attractive system. 
This idea may be understood from the fact that there 
is an ambiguity in calculating the in-medium phase shift 
Eq. ([8]), as it involves a multi-valued ln(a;) function. By 
appropriately choosing different branch cuts, one thus 
may access excited many-body states, in addition to the 
ground state of the system. 

To the best of our knowledge, the proper choice of 
in-medium phase shift was first emphasized by Engel- 
brecht and Randeria in the study of a weakly interact¬ 
ing repulsive Fermi gas in two dimensions in 1992 (^ . 
However, at that time, the connection between attractive 
and repulsive systems was not realized and the concept 
of the upper branch was not established. The meaning 
of the upper branch was only clarified in 2011 by Shenoy 
and Ho, who claimed that by excluding the contribution 
from the paired molecular states in calculating the ther¬ 
modynamics of the system, one could access the upper 
branch of an attractive Fermi gas [s^. This excluded 
molecular pole approximation (EMPA) immediately im¬ 
plies that for the upper branch, the lower boundary of 
the frequency integral in Eq. m should be modified to 
a;s(q), leading to 

1 _ r+oc, 1 

(5H = --^/ yhrep (q,w). (12) 

This expression was later applied by Li and Ho to a 
strongly interacting Bose gas [^. However, despite the 
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clarification of the concept of the upp er branch, in those 
two studies (i.e., Refs. [29| and (58l|h the ambiguity in 
the calculation of the phase shift (5rep(q, w) was not care¬ 
fully treated. The phase shift of the upper branch was 
directly calculated by using 


(q,a; > Ws (q)) = - arctan 


ImT ^ (q, uj) 
ReT”^(Q; w) 


(13) 


without the explanation for the branch cut. Here, the 
function arctan(a;) is the usual inverse tangent func¬ 
tion that takes values in the first and fourth quadrant 
(—7r/2, -l-7r/2) and we have used the superscript “HO” 

to indicate the prescription given by Ho and co-workers. 

It turns out that a more appropriate phase shift for the 
upper branch can be physically defined by the prescrip¬ 
tion 


5rep (q, w) = [5att (q, w) - tt] 0 [w - UJg (q)], (14) 

which can be shown from the viewpoint of the virial ex¬ 
pansion ( 3 ^ . The TT-shift in the above equation can be 
easily understood from the standard scattering theory: 
when a two-body bound state emerges, the two-particle 
phase shift associated with the density of states should 
increase by tt. The prescription Eq. m is therefore 
simply the many-body generalization of the two-particle 
phase shift in the absence of bound states. It should be 
viewed as a physical realization of the EMPA approxi¬ 
mation proposed by Ho and co-workers. 

For a weakly interacting Bose gas (i.e., Og —>■ 0+), The 
two prescriptions for the upper branch phase shift, shown 
in Eqs. (fO)) and ((HI), agree with each other, as a result of 
the large value of ReT”^ (q, w). Towards the strongly in¬ 
teracting limit Og —> -foo, however, the two prescriptions 
differ significantly. In particular, on the BCS side with a 
negative scattering length, the phase shift (5^p(q, w) co¬ 
incides with the phase shift of the ground state branch, 
(5att(q, w). As a result, by changing the scattering length 
and crossing the Feshbach resonance from below, there 
is a sudden branch switch from the upper branch to the 
ground state branch This branch-switching effect 

and the related violation of exact Tan’s relations are 
absent when the more physical prescription Eq. m is 
used. 

In Fig. dKb), we show the in-medium phase shift for the 
upper branch, obtained by performing Eq. (O for the 
attractive phase shift shown in Fig. dja). The Thouless 
criterion ^ (q, w = 0) = 0 is now strictly satisfied. More¬ 
over, for a positive frequency the phase shift becomes 
negative. This leads to a positive fluctuation thermo¬ 
dynamic potential > 0 and a negative pair de nsity 
6n < 0. As analyzed by Engelbrecht and Randeria 
the fact > 0 implies that the ground-state energy 
increases due to the interactions, as it should be for a 
repulsive system. A negative pairing density is also con¬ 
sistent with a repulsive interaction, which, for a specific 
atom, will expel other atoms away from its position, and 
therefore make the effective number density around its 
position smaller. 


In Fig. [Hb), we also report the upper branch phase 
shift at the gas parameter na^ = 0.01 (grey crosses) and 
nOg = — 1 (green stars). It is worth noting that the neg¬ 
ative value of the gas parameter (i.e., on the BCS side 
above the Feshbach resonance) actually means stronger 
repulsions between atoms, as indicated by the large ab¬ 
solute value of the phase shift. In contrast, for a posi¬ 
tive gas parameter, the interaction effect becomes weaker 
with decreasing the gas parameter. 


C. Large-A expansion 


The generalized NSR approach was used earlier by Li 
and Ho to investigate a strongly interacting Bose gas 
near unitarity [^ . A large mechanically unstable area 
was found when the temperature of the system is be¬ 
low T < 5Tco 'iTp, which renders the approach useful 
only at extremely high temperatures. Here, we show that 
the mechanical instability is artificial and caused by the 
inappropriate treatment for the strong pair/density fluc¬ 
tuations in the NSR approach. It can be cured by the 
so-called large-iV expansion technique 

In the large-A expansion, we assign an additional fla¬ 
vor degree of freedom to bosonic atoms (z, j = 1, • • • , N) 
and thereby extend the model action to, 


5 = 


dTd:x. 


. 1=1 ^ 



+ 


2N 


N 


ij=l 


(15) 


By introducing a pairing field 
U ^ 

(x, r) = ^ ^ (x, r) z/ji (x, r) (16) 

i=l 

and again decoupling the interatomic interaction via the 
standard Hubbard-Stratonovich transformation, we inte¬ 
grate out the atomic fields and obtain the grand thermo¬ 
dynamic potential per flavor 

§ = f!„ + lzn + o(T). (17) 

up to the first non-trivial order of 0(l/iV) Here, 

for the metastable upper branch, Hq and <517 are given by 
Eqs. (0 and ([12]), respectively. 

It is clear that in the large-A expansion we have intro¬ 
duced an artificial small parameter 1/A, which can be 
used to control the accuracy of the theory of strongly in¬ 
teracting Bose gases. The NSR approach, which is based 
on the summation of infinite ladder diagrams [HJ, 
should be understood as an approximate theory obtained 
by directly setting N = 1. However, such a procedure 
cannot be justified a priori in the strongly interacting 
regime, as the controllable parameter 1/A is already at 
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FIG. 2: (Color online) Energy as a function of the artificial 
controlling parameter l/N at T = 2Tco (black squares for 
nal = 10~® and red circles for na^ = 1) and T — lOTco 
(blue crosses with na® = 1). The energy is measured in units 
of the energy Eq of a non-interacting Bose gas at the same 
temperature. The lines are linear fits to the small l/N data. 
Note that, at low temperatures and strong interactions (i.e., 
red circles), we are not able to find solutions for N — 1 because 
of strong correlations. 


0.039 0.181 0.840 ka 

F s 



the order of unity. Indeed, the appearance of the large 
mechanically unstable regime at low temperatures, men¬ 
tioned at the beginning of this subsection, is precisely 
an indication of the breakdown of the procedure of di¬ 
rectly setting N = \. A more reasonable treatment is 
to first solve the thermodynamics of a A^-flavor system 
with iV ^ 1 and then linearly extrapolate all the desired 
physical quantities - as a function of 1/iV - to the limit 
of = 1. This large-expansion idea has been suc¬ 
cessfully applied to a strongly interacting two-component 
Fermi gas in the unitary limit (dol . |4l| . The equation of 
state and the Tan contact near the quantum critical point 
/i = 0 was then accurately predicted [d^. In this work, 
we anticipate that the same large-A^ expansion technique 
could also lead to very useful information for a unitary 
Bose gas in the quantum critical region. 

In Fig. [21 we show the 1/A^-dependence of the total 
energy of an interacting Bose gas at different gas param¬ 
eters and temperatures, obtained by solving the coupled 
equations Eqs. and CZD, and subject to the 

number equation h/N = n = uq + Sn/N for the num¬ 
ber density per flavor n/N. At weak interactions (black 
squares) or high temperatures (blue crosses), roughly the 
energy changes linearly as a function of 1/A^. The linear 
extrapolation approximation used in the large-A^ expan¬ 
sion therefore does not make signihcant difference. How¬ 
ever, for a strongly interacting Bose gas at relatively low 
temperatures (red circles), the dependence is highly non¬ 
linear. In particular, we are not able to find physical 
solutions when the number of flavors N < 2. Therefore, 
it becomes crucial to keep only the linear term in the 
1 /A^ expansion, which provides the first non-trivial and 


FIG. 3: (Golor online) The energy (a) and pressure (b) as 
a function of the gas parameter na® at T = 2Tco (black 
solid lines) and T — 4Tco (red dashed lines), normalized re¬ 
spectively by their corresponding results of an ideal, non¬ 
interacting Bose gas at the same temperature. The results 
from a path-integral Monte-Carlo calculations are also shown 
1^ . with squares for hard-sphere potential and circles for 
soft-sphere potential. 


non-pertrubative knowledge about a strongly correlated 
many-body state. 


III. RESULTS AND DISCUSSIONS 

In this section, we present our large-A^ results, cal¬ 
culated by the linear extrapolation towards the limit 
1/A^ = 1. In practice, we solve the generalized NSR 
approach with A^ = 50 — 100 for the chemical poten¬ 
tial and the energy E{N), and then expand them 

around the corresponding non-interacting values po and 

Eo, 

^i{N) = ^^o + S^^/N + 0{l/N^), (18) 

E{N) = Eo + SE/N + 0 (l/N^), (19) 

to extract the corrections Sfi and 6E. This leads to the 
large-A^ expansion results ^ = ^o + 6fj, and E = Eq + SE. 
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FIG. 4: (Color online) Two-body contact I 2 as a function of 
the gas parameter na^ at T = 2Tco (black solid line) and T = 
dTco (red dashed line). The result from the zero-temperature 
Bogoliubov theory is shown by the blue dot-dashed line. 


A. Crossover to strong repulsions 
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In Figs. m^) and ISKb), we present respectively the 
energy and pressure of an interacting Bose gas at two 
temperatures T = 2Tco (black solid lines) and T = 4Tco 
(red dashed lines), as a function of the gas parameter na^, 
or kpUs if we convert the number density n to a Fermi 
wavevector kp = (Stt^u)^/^. The large-7V expansion re¬ 
sults are compared with available path-integral Monte 
Carlo calculations for a hard-sphere (squares) and soft- 
sphere potential (circles) jH]. For weak interactions (i.e., 
naf = 10 “® and 10 “^ or kpag < 0 . 2 ), our predictions 
agree well with the ab-initio simulations. For strong in¬ 
teractions with strength kpOg 0 . 8 , there is a significant 
difference. This is due to the effect of non-negligible ef¬ 
fective range of interactions tq used in the Monte Carlo 
simulations (i.e., |A:f?’o| l)i which leads to a sizable cor¬ 
rection to the energy and pressure. In our calculations 
with a contact interaction, the range of interactions is 
strictly zero. 

In Fig. m we show the evolution of the corresponding 
Tan’s contact with increasing gas parameter. Tan’s con¬ 
tact measures the density of pairs at short distance and 
determines the exact large-momentum or high-frequency 
behavior of various physical observables [l^. It therefore 
serves as an important quantity to characterize a strongly 
interacting many-body system. In particular, experimen¬ 
tally it can be measured from the momentum distribution 
[ 12 , [ 1 ^, which takes a k~'^ tail in the short-wavelength 
limit, i.e., n(k) —>■ I/k^. At hnite temperatures, the con¬ 
tact can be conveniently calculated by using the adiabatic 
relation 


T2 = 


~w 


do. 


da. 


-1 


■iTp 


( 20 ) 



r/ro 


FIG. 5: (Golor online) Temperature dependence of the chem¬ 
ical potential (a), energy (b) and entropy (c) of a unitary 
Bose gas. For comparison, we show the second-order virial 
expansion predictions by red empty circles and the ideal gas 
results by dot-dashed lines. The latest QMG results at zero 
temperature are also plotted by using stars in brown [35l| . 


We have used the subscript “2” to emphasize the fact that 
in our calculations we do not consider the three-body Efi¬ 
mov physics and the associated inelastic collisions. These 
effects are instead captured by a three-body contact I 3 , 
which can be defined through an adiabatic relation for 
a three-body parameter. We refer to Ref. [s^l for more 
detailed discussions. 

The two-body contact is an increasing function of the 
interaction strength. At small gas parameters, our results 
are in good agreement with the weak-coupling predictions 
of a zero-temperature Bogoliubov theory (thin blue dot- 
















dashed line) [^ . 
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21bog = (dTrnas 



( 21 ) 


The slight increase in our large-TV expansion results is due 
to the finite temperature effect. At large gas parameters 
nOg, the contact tends to saturate to a universal value 
that depends only on the temperature, as it should be. 


B. Unitary Bose gases 


We are now in position to discuss the universal ther¬ 
modynamics of a unitary Bose gas. In Fig. [5l we present 
the chemical potential, energy and entropy, as a func¬ 
tion of temperature. For comparison, we also plot in 
dot-dashed lines the temperature dependence of an ideal, 
non-interacting Bose gas. For the chemical potential 
and energy, our results lie systematically above the non¬ 
interacting results, clearly indicating the consequence of 
strong repulsions. They tends to converge to the zero 
temperature quantum Monte Carlo predictions (brown 
stars) with decreasing the temperature (35| . In contrast, 
the entropy seems to be less affected by strong repulsions. 
The insensitivity of entropy on the interatomic interac¬ 
tions was also previously found for a unitary Fermi gas 


At high temperatures with a small fugacity z = <C 
1 , we may use the virial expansion theory to study the 
universal thermodynamics |48| . For a unitary Bose gas, 
the virial expansion of the grand thermodynamic poten¬ 
tial takes the form. 


n = no-^{z^Ab 2 + z^Ab 3 + ---), (22) 



TIT ^ 

cO 


FIG. 6: (Color online) Temperature dependence of the two- 
body contact X 2 in the unitary limit. The prediction from 
a second-order virial expansion, Eq. ipsii. is shown by red 
empty circles. At zero temperature, the brown star and green 
solid circle refer to the latest QMC result (ssll and the number 
extracted from the recent measurement at JILA for a trapped 
unitary Bose gas [^. Is^l . The blue cross indicates the zero- 
temperature contact of a unitary Fer mi g as, measured very 
precisely by using Bragg spectroscopy [6^ . 


a direct consequence of the adiabatic relation, we have 
the expansion. 


^ Sirm ksT , „ 

I 2 = — {z C 2 + z C 3 • • •) , (24) 


h? X 


dB 


where c„ = \^g{dAbn/daJ^) is the so-called contact co¬ 
efficient [ 6 l|. For a unitary Bose gas, by using Beth- 
Uhlenbeck formalism it is easy to show that C 2 = 2 / 7 r. 
Using Eq. (^51) . we then find 


where XdB = [27r/i^/(mfesT)]^/^ is the thermal de-Broglie 
wavelength, and A &2 = — is the second-order virial 
coefficient for strong repulsions [d^, H^, which can be 
easily calculated by using Beth-Uhlenbeck formalism . 
Up to the second order, we can solve Eq. together 
with the number equation Eq. (nnj. Eor the fugacity, we 
find that 


lie 

V ^ / 


(23) 


at T ^ Tp. The virial predictions for the equation of 
state are shown in Eig. [5] by red circles and agree well 
with the large-A expansion results at high temperatures. 

In the unitary limit, we may calculate the universal 
contact by using the adiabatic relation, Eq. (j20L shown 
in Eig. [HI With increasing temperature, the contact ini¬ 
tially increases and then decreases, giving rise to a peak 
structure at the temperature T ^ dT^o. The decrease of 
the contact at high temperatures can be well understood 
by using the virial expansion theory for the contact. As 


h _ 

nkp 3 \ T J 


(25) 


Thus, at high temperatures the contact decreases as T~^. 
There is no apparent physical explanation for the increase 
of the contact at low temperatures. However, we notice 
that with decreasing temperature toward zero tempera¬ 
ture, our large-A expansion result seems to be consistent 
with the zero-temperature value predicted by the latest 
quantum Monte Carlo simulation [s^. 

Eor comparison, we show the experimental data of the 
contact I 22 I I , a nalyzed by Smith and co-workers (green 
solid circle) |32l | . The significant discrepancy between ex¬ 
periment and our large-A theory should be largely due to 
the unknown temperature in the experiment, as the Bose 
cloud could be significantly heated by atom losses (^ . 
We also show the zero-temperature contact of a unitary 
Eermi gas (blue cross), which has been both calculated 
and measured very accurately (^ . It is interesting that 
both the unitary Bose and Eermi gases have similar con¬ 
tact at zero temperature, indicating that a 3D Bose gas 









may also have the tendency of being fermionized at strong 
repulsions, analogous to a Bose gas in one dimension. 

IV. CONCLUSIONS 

In summary, based on the upper branch idea and large- 
N expansion technique, we have developed a unified the¬ 
ory for a normal, strongly interacting Bose gas. The 
theory reproduces the path-integral Monte Carlo sim¬ 
ulation in the weak-coupling limit [s^. While at high 
temperatures, it nicely recovers the known results from 
a quantum virial expansion calculation (^ . |48| . Thus, 
we anticipate that the universal thermodynamics pre¬ 
dicted by our theory could be qualitatively reliable. A 
useful check may be provided by experimentally measur¬ 
ing the finite-temperature contact of a unitary Bose gas 
through momentum distribution or momentum-resolved 
radio-frequency spectroscopy ■ 


Our results complement the earlier studies of a con¬ 
densed strongly interacting Bose gas. It is worth nothing 
that our theoretical framework can naturally be extended 
to include the condensation (i.e., T < Tco) by using a gen¬ 
eralized Nozieres-Schmitt-Rink approach below the su¬ 
perfluid transition temperature |56l |. This extension will 
be addressed in a future investigation. 
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